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1 Introduction 

Let the deterministic motion of a passive particle in a random velocity field V(t,a;) be 
described by 

^4v(^,..W) (1, 

where e is a small parameter and V(t, x) = {Vi{t, x),- ■ ■ , Vd{t, x)) a time-stationary, space- 
homogeneous, incompressible, V ■ V = 0, zero-mean velocity field. The question is to deter- 
mine the asymptotics, as e tends to zero, of the displacement x^{t) from the statistics of the 
velocity. 

Scaling like (0) arises naturally when the velocity has a non-zero mean drift and a small 
fluctuation V + e\'{t, x). By changing to the moving coordinate system x ^ x -\- Vt and 
taking the long time limit t -^ t/e^ we have a equation of the form (0) with V(t, x) = 
V'(t,x + Vt). 

Under the assumption of strong mixing in time among other things, Kesten and Papan- 
icolaou proved that Xi;{t) converges to a Brownian motion with the diffusion coefficient 
given by a Taylor-Kubo formula 

oo 

D*^ = J {E[V,it, 0)V,(0, 0)] + E[V,{t, 0)V,(0, 0)]} dt (2) 



(see also Carmona and Fouque p for the corresponding result of stochastic flows). No 
assumption on the decay of velocity decorrelation in space is required. 

Convergence to Brownian motion in the absence of molecular diffusion is often referred 
to as turbulent diffusion because the large scale diffusive motion is a result of random flows 
instead of molecular diffusion. 

Does the turbulent diffusion theorem hold for velocities lacking the temporal mixing prop- 
erty but with correlations decaying in space? To study the interplay of temporal and spatial 
correlations in a precise way, we consider the class of time-stationary, space-homogeneous, 
isotropic, Markov, Gaussian velocity flelds V(t, x), deflned on the probability space (f2, V, P), 
with two-point correlation tensor R = [Rij{t, x)] given by the Fourier transform 

R,,{t,x) = E[V{t,x)0V{O,O)]= (3) 

/ cos (k ■ x)e-l''l"'*R(k)dk 

where E stands for the expectation. The spectral density R(k) is given by a power law 

ftW-n^^fl-^). "<1. .5>0. (4) 



with a compactly supported, continuous function a : [0, +00) -^ R^. The factor I — ^^ in 
(I) ensures that the velocity field is incompressible. 

The function exp (— |kp^t) in (|^) is called the time correlation function of the velocity V. 
The spectral density R(k) is integrable over k for a < 1 and, thus, (^-(^ defines a velocity 
field with a finite second moment. The ultraviolet cut-off K is needed to avoid divergence of 
integral over large |k|. Because of the ultraviolet cutoff V is jointly continuous in t,x and 
is C°° in X almost surely. 

The parameter a is directly related to the decay exponent of R : R(ic) ~ |a3|"^^ for 
|£c| ^ 1. As a increases to one, the spatial decay exponent of R decreases to zero and, 
consequently, spatial correlation of velocity increses. On the other hand, for j3 > 0, the 
velocity field lacks the spectral gap and thus strong mixing property (see [0). We restrict 
our attention to the case a < 1, (3 > which corresponds to a velocity field with arbitrarily 
long scales but not the strong mixing property. 

What is the region in the (a,/3) plane where the turbulent diffusion theorem, with the 

Taylor-Kubo formula, holds? It is easy to find the necessary condition by imposing the 

convergence of the Taylor-Kubo formula. A plain calculation 

00 
D;, = fR,,{t,0)dt (5) 

I {5., - ^)ai\k\)\k\'-''^~'' I expi-\k\'f't)dtdk (6) 

Rd 



ij J "^*J 





R'i 



a fc^fcj x Q(|k|) dk 



leads to the condition 



a + p <1. 



As it turns out, (^ is also sufficient. This is our main result which generalizes the turbulent 
diffusion theorem of |0] to a class of non-mixing flows. 

Theorem 1 Let V(t, x) be a Markov, Gaussian velocity field with correlation given by (Qj- 
(^. Then, for a + f3 < 1, the invariance principle holds for the displacement x^{t). Namely, 
x^{t), as continuous process, converges weakly, as e I 0, to a Brownian Motion with covari- 
ance matrix given by formula (j^). 

Additional molecular diffusion in (|l|) would act like a regular perturbation to the eddy 
diffusivity defined by the Taylor-Kubo formula. 



What happens if (^) is violated? From the divergence of the Taylor-Kubo formula one 
judges that the long time asymptotics should be faster than diffusion and, consequently, a 
different scaling limit is necessary. It turns out that under an anomalous scaling the dis- 
placement converges to a fractional Brownian motion. The fractional Brownian motion limit 
theorem will be studied in a forthcoming paper. The dichotomy is conveniently represented 
in the figure. 

2 Multiple stochastic integrals 

By the Spectral Theorem (see, e.g., [|l^) we assume without loss of generality that there 
exist two independent, identically distributed, real vector valued, Gaussian spectral measures 
Vi(t,-), / = 0,1 such that 

V{t,x)= fVo{t,x,dk), (9) 



where 

Vo(t, X, dk) := co(k ■ x)Vo{t, dk) + ci(k • cc)Vi(t, dk) 

with co(0) = cos (0), ci(0) = sin (0). Define also 

Vi(t, X, dk) := -ci(k ■ a:;)Vo(t, dk) + co(k ■ x)Vi{t, dk). 

We have the relation 

d\'o{t,x,dk)/dxj = kjVi{t,x,dk), (10) 

dV^{t,x,dk)/dxj = -kj\'oit,x,dk). (11) 

One can check that / Vi(t, a;, dk) is a random field distributed identically to and indepen- 
dently of V. 

We define the multiple stochastic integral 

J ■■■ J ip{ki,---,kf^)Vi^{ti,Xi,dki) 'S)---<»Vij^itN,XN,dkN) (12) 

for any li,- ■ ■ ,lj\i G {0, 1} and a suitable family of functions ip by using the Fuhini theorem 
(see (|l^) below). For t/^i, ■ ■ ■ , ipjq G S{R'^), the Schwartz space, and /i, ■ ■ ■ , /a? G {0, 1} we set 

/ ■•■ / ipi{ki)---ipN(kN)yi^iti,Xi,dki) ®---'S)Vi^{tN,XN,dkN) (13) 

:= I 'ilJi{ki)\'i^{ti,xi,dki) (g)---(g) / V7v(k7v)Vij^(tAr,ccAr,rfkiv). 



We then extend the definition of multiple integration to the closure 7i of the Schwartz space 
S{{R'^)'^ ,R) under the norm 



E 



Vi^{ti,xi,dki] 



:=|---|V^(ki,---,MV^(k;,---,k^) 
Yij^(tN,XN,dkN) ■■ Yij^(ti,xi,dk[) (g) • • 



(14) 



Yij^(tN,XN,dk'j 



Ni 



The expectation is to be calculated by the formal rule 



E 



Vi^i{t, X, dk)Vv^i\t\ x\ dk!)] = e-l''l"'l*-*'l5;,//Co(k ■ [x - x'))Ri^i,{k)5{k - k!)dkdk! . 



Among various possible definitions of multiple stochastic integral (see, e.g., [|r^, 13, I^|) this 



approach (|jT5[) seems to best serve our purposes. 

We denote the stochastic integral ([T2|) by '^iiti, ■ ■ ■ , t^v, Xi, • ■ • , X]^) or, sometimes, simply 
by ^i. Here I = (/i, ■ ■ ■ , /at), /i, ■ ■ ■ , /^v ^ {0)1}- Sometimes when i = {ii,- ■ ■ ,id), with 
ill' ' ' I'i'd ^ {1) 2, ■ ■ ■ , (i} is fixed we shall denote the corresponding component of the above 
tensor valued random field by ^i^. 

Note that ^E^^^i G if^(V) - the Hilbert space obtained as a completion of the space of 
A^-th degree polynomials in variables / '?/'(k)V(t, a;, k) with respect to the standard L^ norm. 



Proposition 1 For any (ti,Xi),- ■ ■ ,{tiy,xiy) E R x W^ and p > 0, "^li belongs to L'P{n) 
and 

(E|vl>,,in^/^<C(E|M/,,ip)'^' (15) 

with the constant C depending only on p, N and the dimension d. Moreover, "^^i is differ- 
entiable in the mean square sense with 



Vxj'^e,i{ti,---,tN,xi,---,XN) = J --' J kjilj{ki,---,kN) 
Vi^,i^{ti, xi, dki) ■ ■ ■ Vi_i^^i^{tj, Xj, dkj) ■ ■ ■ t^^,i^(tAr, Xn, dk^) 



(16) 



The proof of Proposition ^ is standard and follows directly from the well known hyper- 
contractivity property for Gaussian measures (see, e.g., 0, Theorem 5.1. and its corollaries), 
so we do not repeat it here. 

The field V is Markovian 



E 



tpik)Vi{t,x,dk)\ V_ 



-lkp/3(t-s) 



^{k)yiis,x,dk), 1 = 0,1, (17) 



for all ip G S{R'^, R), where Va,b denotes the a-algebra generated by random variables V(t, x), 
for t G [a, b] and x G R'^. 



To calculate expectation of multiple product of Gaussian random variables, it is con- 
venient to use a graphical representation, borrowed from quantum field theory. We refer 
to, e.g., Glimm and Jaffe |Q and Janson 0]. A Feynman diagram T (of order n > and 
rank r > 0) is a graph consisting of a set B[T^ of n vertices and a set E{J^) of r edges 
without common endpoints. So there are r pairs of vertices, each joined by an edge, and 
n — 2r unpaired vertices, called free vertices. B{J^) is a set of positive integers. An edge 
whose endpoints are m,n & B is represented by ffm (unless otherwise specified, we always 
assume m < n); and an edge includes its endpoints. A diagram JF is said to be based on 
B{T). Denote the set of free vertices (i.e. points which are not endpoints of any edges) by 
A(^), so A{J^) = JF \ E{J-'). The diagram is complete if A{J-') is empty and incom,plete, 
otherwise. Denote by G{B) the set of all diagrams based on B, by Qc{B) the set of all 
complete diagrams based on B and by Qi{B) the set of all incomplete diagrams based on B. 
A special class of diagrams, denoted by Qs{B), plays an important role in the subsequent 
analysis: a diagram JF of order n belongs to Qs{B) if Ak{J^) is not empty for all A; = 0, 1, ..., n. 
Let B = {1, 2, 3, ..., n}. Denote by J-'\k the sub-diagram of JF, based on {1, • • • , k}. Define 
Afc(^) = A(^|,). 

We work out the conditional expectation for multiple spectral integrals using the Markov 
property (|T^). 

Proposition 2 For any function if) ^Ti and /i, ■ ■ ■ , /at G {0, 1}, ii,- ■ ■ ,iN G {1, ■ ■ ■ , d}, 



E 



tpiki, ■ ■ ■ ,kN)Vi^,i^{t,Xi,dki) ■ ■■Vi^^i^{t,XN,dkisi) \ V. 



oo,s 



:i8) 



Y, /■■•/exp<- Yl Ikmp'^(t-s) >V^(ki,---,kAr)Vs,Xi,-,a3Ar(c?ki,---,(ikAr;JF) 

with 

Vs,x^,...,Xi^{dki,---,dkN;J^) := H ^im,im{s,Xm,dkm) (19) 

Y[ \^^ _ ^-i\^^\-^M^.n(t-s)-]^ E[Vl^,Js,Xr,,dkr.)Vl„,^{s,Xr„dK)]. 

Proof. Without loss of generality we consider '?/'(ki, ■ ■ ■ , k^v) = l^^(ki) • ■ ■ lAj^(k]\f) for 
some Borel sets Ai,- ■ ■ , Ajsf. 

Vi{t, Ai) = Vf (t, Ai) + V/(t, Ai) where V,°(t, ■), V;^(t, ■) are the orthogonal projection of 
Vi(t, ■) on the space L'^^o s ^^^ i^^ complement, respectively. Here L^^ denotes L^ closure of 
the linear span over V(s, x), a < s < b, x & R^. The conditional expectation in (pIS]) equals 



E n WLlSt^^rn)Vl,St^Ar.)\ \{ \/,l, Jt, A 



6 



The statement follows upon the application of the relations 



V°(t,A)= I e 

A 



and 



E 



\l{t,A)0ylit,B) 



J J 6i^v {e [Vi(t, rfk) ® Vp(t, rfk')] - E [Vf (t, c/k) ® V°(t, (ik')] } . 



A B 



For fixed x ^ R^ we define a process 

^£(t; as) := ^<.(t, ■ ■ ■ , t, cc, ■ ■ ■ , a;), for t > 0. 

Let us introduce, after Chapter 3 of p, the pseudogenerator oi the above process C(t;x) 
t > as 



lim 

510 



(5 



in the L^ sense. After elementary calculations using ([18|) we obtain that 

C{t; X)=Y.J--- y'(|kp|2/5+|m2/^)V^(ki, ■ ■ ■ , \iN)KK-K)K,k^M^^^ a;, dki, • • • , 0?kAr)rfkprfk, 

(20) 
Here the summation extends over all edges pq with vertices belonging to the set S = 
{1,---,A^} and Vri(s, a;, rfki, ■ ■ ■ , t/k^v) indicates that the multiple stochastic integration 
is taken over the product of stochastic measures V;^(s, x, (ik„) for all n G 5 \ {p, q}. 
The fundamental property of C{t; x), t > is that for any fixed x 



M{t] x) := ^^(t; x)- j C{s; x)ds 



(21) 



is a continuous trajectory martingale, cf. 0. In fact we can give an explicit expression on 
its quadratic variation process < M >t, t > in terms of the spectral representation. 



Proposition 3 For any fixed Xi,- ■ ■ , Xjq G R^ 



< M(,-;a;) ><= fw{s,x)ds 



(22) 



where 



W^(t,a;):=2X:/---//---/^(ki,---,k^)V'(k;,---,k'^)a(|k,|)|k,|2/^+2-2"-'^ 



pq 



( I j, p^ J O V— (t, X, dki, ■ ■ ■ , dk2Ar)5(kp - kq)6i^^i^dkpdkqdkpdkq. 

Here we adopt the convention that k^ = kN+n, In = ^N+n- The summation extends over all 
edges pq with vertices p, q belonging to the sets {1, ■ ■ ■ , A^} and {N + 1, ■ ■ ■ , 2N} respectively. 

Proof. The quadratic variation of M(-; x) coincides with the pseudogenerator of M(-; x) ® 
M(-; x). After an elementary calculation we get that the left hand side of ( p2D equals 

C'{t; x) - -^e{t; x) (g) C{t; x) - C{t; x) ® *^(t; x), 

with C'{t; x), t > denoting the pseudogenerator of \E'£(t; x) ® ^^(t; x), t > 0. The proof of 
the proposition can be then concluded from an application of formula (^) I 

3 Proof of Theorem J 

We begin this section by defining the key concepts: the n-th X-corrector, X\\ ^-^id the n-th 
X-convector, V^ recursively: 

(23) 
ds (24) 

(25) 

The differentiation in (^) is taken in the mean square sense. 

The definitions requires the integrability of convectors and the differentiability of correc- 
tors, as stated in the following lemma. 

Lemma 1 Xx > '^x ^^^^ L^ integrable for any p > 1, X > and n > and satisfy the 
following hypercontractivity property 

{E\xt\t,x)\f' <C{E\xt\t,x)\^y^' (26) 



VTis,x) - 


= Y{s,x) 






+ 00 


X^x\t,x) :- 




vt'\s,x) v_oo,; 


U?Ht,a.) :- 


= Vis,x)-Wxi 


''\t,x). 



{E\lj'^-\t,x)\^)''<C{E\lJ^-\t,x) 



where the constant C > depends on n,p, d but not on X. 

There exist versions of Xx ^'^'^ ^x 
continuous in both (t,x), C°° smooth in x 



There exist versions of Xx ^'^'^ ^x "which are jointly stationary in the strict sense, 



L^ integrability follows from the hypercontractivity property on Gaussian Measure spaces 
(see Theorem 5.1, Nelson's theorem and the corollaries in 0). The existence of regular 
versions of the fields can be proved by using aforementioned Nelson's theorem in conjunction 
with Kolmogorov's criterion of continuous processes (see Theorem 3.2.5 and corollary, |]I|]). 

This and other lemmas in this section are proved in the next section. 

We define the n-th rescaled corrector along the path by 

Xi")(t):=5"x5^(^,^.(t)) (27) 



and the n-th rescaled convector by 

t 

Because of the incompressibility of V, both xi"H^) ^-^d U^"'-*(t) are stationary (see Theorem 



Uf(t):=e"-iu5)(^,a;,(t)). (28) 



2of 0). 

Next we represent Xx by multiple spectral integrals. To simplify the presentation we 
introduce some notations. For m = {mi, ■ ■ ■ , nip) and p' < p we set m|p/ := {mi, ■ ■ ■ , nip'), 
mlP := {nipi, ■ ■ ■ , nip) and (m, k) := {mi, ■ ■ ■ , nip, k). 

We introduce next the proper functions of order n,(T : {!,■■■ ,n} ^ {0, 1} that appears 
in the statement of the next lemma. The proper function of order 1 is unique and is given 
by cr(l) = 0. Any proper function, cr', of order n + 1 is generated from a proper function cr 
of order n as follows. For some p < n, 

(29) 



(T'{n + 1) 


:= 


a'{k) 


:= (T{k) for k < n and k ^ p 


^'{P) 


■■= i-t(p). 



(30) 

In other words, each proper function cr of order n generates n different proper functions of 
order n + 1. Thus, the total number of proper functions of order n is {n — 2). 

Lemma 2 For any ii G {!,■ ■ ■ ,d} and X > we have 

xtlit.x) (31) 

= E (■■■Ut^{\^i.---.K)Pn-i{:F) y\ n K„.Jt,a;,rfM 

where Lp{)j are some functions with sup \^icr\ < 1 ciiT'd 

V Ik I 

n-l ^ I '"I 

9 



Q^^) = _n E[K„,„(o,rfk^)v;^„,^,(o,rfk^O] (33) 

where cr is summed over all proper functions defined over {1, ■ ■ ■ , n}. 

Remark. Indexing terms of the sum in (|3l|) with proper functions keeps track of differ- 
entiations we have to perform in the process of constructing consequtive ^'^"'^-s. The idea is 
as follows. The spectral representation of the n-th order corrector contains terms which are 
multiple spectral integrals with respect to the product measures of the form 

n K,^,iJt,x,dkJ. (34) 

Sigmas stand for the fact that there are two types of measures Vq, Vi involved which are 
orthogonal to each other. In order to construct the next generation corrector we have to 
apply nabla operator in x to ( P^ using the Leibnitz rule. This procedure consists in picking 
some factor Vap,ip(t,x,d'kp) for a certain p < n, belonging to An^J-") and changing that 
term to Vi-a ,« {t,x,d\<.p). The factor kp arising in the process will be absorbed into the 
remaining expression. Finally to form ;i^("+^) we have to multiply yet the entire expression 
by VoA„+i(t, X, dkn+i). The manipulations we have just described are reflected by the rule of 
forming the proper function of the consequtive generation described by formula (p9|). 
The utility of the correctors and convectors depends on the following lemma. 

Lemma 3 We have the recursive scheme, for n > 1, 

t t t 

jvf~'\s)ds = -xfit) + xf (0) + Jx'e'\s)ds + Ju^:^\s)ds + Mi^-'\t). (35) 



where M^^^ (t) is a Brownian Motion with the covariance matrix 

f |kpM|k|) f, k®k\ 

D. = J |kp...-.(|i,|2/3 + ,.). (I - -J^ dk. (36) 

and M^^\t),n > 1, are continuous martingales with respect to the filtration V_ 



^'77 



For a + f] < 1 there exists '~f{n) > such that 

lim£-^(")E|xi"H^)P = 0, forn>l (37) 

eiO 

lim£-'^(")E|U("H^)P = 0, for n > 2[(3] + 1, [(3] : the integer part of (3 (38) 

eiO 

limE|£-^(")Mi")(t)|P = 0, forn>l (39) 

eiO 



10 



We apply now Lemma 3 to prove Theorem |T]. By (|3^) we have 

^,(t) = Mf(t) + /?,(t) (40) 

with the remainder Relt) 

Re{t) = i?,,,(0) - Rs^eit) + Rrn,e{t) + RaA^) (41) 

where 

2[/3]+l 
n=l 

is stationary, 

RaAt)--= E j xf\s)ds + j Vf^^^+'\s)ds (42) 

"=1 

additive and 

2[/3] 

i2„,,(t):=^Mi")(t) 

n=l 

a vanishing martingale. 

Clearly, M}°''(t) converge, as continuous processes, to the Brownian motion with the 
covariance Dq = lime^oDg (see, e.g., [^). 

To finish the proof it remains to show tightness limP[ sup |-Re(t)| > (^] = 0, VT > 

eiO 0<i<T 

0, 5 > 0. By Lemma |^ it suffices to show limP[ sup \Rs£{t)\ > 5] = 0. Since Rgs is 

£-10 o<t<r 
stationary, 

P[sup |i?,,,(t)|>(5]<^P[ sup |P.,e(t)|>(5]. (43) 

0<t<T S 0<t<Te^ 

Using (|40|), (|4l|) we bound the right side of (^) by 

^P[ sup [\RaAt)\ + I^.WI] > <^] + ^P[ sup [|P™,.(t)| + |Mf (t)|] > 6]+ (44) 

^ 0<t<Te2 e 0<t<Te^ 

|P[|P.,,(0)|><5]. 
The last term in (|4^ ) is smaller than 

E|P.,,(0)|PT 

by Chebychev's inequality. 

Let p be such that 7p > 4 where 7 = min 7(n). Expression (H3) tends to as e ], by 

l<n<A'' 
(0)- 



(45) 
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To prove that the second term of (44) vanishes as e | we only need to show that for 
any n > 1 



T 



P[ sup |Mf)(t)|] >5] 



£ 0<t<Te'^ 

tends to as £ I 0. Using (^) and Kolmogorov's inequahty we bound (|46|) by 



(46) 



E|MJ")(T)|PT 

which vanishes with e. 

The first term of (|^) can be bounded by 

'^2rp 2[/3] + l 



<C 



%n{n)~2 



TdP 



C 



S^ 



E'E|x["H0)r + ^E|.^[/^]u(f]+^)(0,O)r+^P[/|V(s,.=r(.))|rf.>^ 

n=l " ^ n ^ 



(47) 



Expression (^) tends to zero by Lemma |^ and that 



hm sup E 

£j.O 



T -iP 

'V(s,ex(s))\ ds 



LO 



1 

< TP-^ fE\V{s,ex{s))fds 



= E|V(0,0)|^TP<+oo. 
For p > 2 the last term in (^Tf ) vanishes with e. Thus the proof of Theorem |l] is complete. 

4 Proof of Lemmas 2 and 3 

Proof of Lemma |2|. We prove (|3l|) by induction. The case n = 1 is obvious. We choose 



,(0) 



and 



iflj = 1. Suppose that the result holds for n. Thus by (| 

+ 00 

= J e-^(^-*)E {V(s, x) ■ VxS {s, x) I V_oo4 c^s 



I -l-cxj 

E E //■•■/e-^(-V[3-(kil,---,k„)P„_i(-F)g(^) 



1 






Vo(s, a;, rfk„+i) ■ V <( H ^fT^.i^l^, a;, rfk^) ^ c^s|V- 



oo,t 



(48) 



To calculate (H) we decompose each t4-^j(s, as, dk) as t4-,i(s, a;, dk) = t^j(s, a;, (ik) + 
^'^^(s, a;, (ik) where 

V;°,(5, X, cik) = e-l''l'''(-*)K,(t, a;, dk) (49) 
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is the orthogonal projection onto L?i^ ^ and 

V^],{s, X, dk) = (1 - e-l'^l''(^-*))K,(t, X, dk) 

the orthogonal complement. We have Vo-(s, x, dk) = V°(s, x, dk) + V^(s, x, dk). 
Expression (^81) becomes 



(50) 



-too 

Y: / E [I ■ ■ ■ I e-^(^-*V[3-(ki, ■ ■ ■ , k„)P„-i(^)Q(^)/C(^) I V_oo,t 



ii2,:?^,cr t 



ds (51) 



with 
/C(^) := 



1 



X+ T \kW ^ Vo^"-(.,x,dk„+0-V<| n V^':,Js,x,dK 

Z^ I nil e={e,} 



,mGA„(:r) 



The term corresponding to Qj = 1 vanishes, as is shown in the following calculation 



E 



/ /■■■/e-'^'"*Vl3-(ki,---,k„)P„_i(^)Q(J-) A+ E Ik^l^'^l (52) 



ds 



Vlis, X, c/k„+i) ■ V n Kl^Js, X, dk„ 

(+00 
E / [/■■■/e"^(-*)¥.S3.(ki,---,k„)P„_i(^)g(^) 



-1 



A+ ^ \km\'n Yi{s,x,dK^,) n yLijs,x,dK 

= 

by homogeneity. 

We use eqs. (|IU[) and ( [TT| ) to compute 



ds 



yo{s,x,dK+i)-V I n '^^m,i^is,x,dkrr 

^ k„/ ■ V w (s, a;, rfk„+i) ]^ t/^„,' , (s, a;, dk„ 



-^n+1 



(53) 



where cr"* is a proper function obtain from a via procedure described by ( pOD with p = m'. 
By (H), dg), (H) and the definition (|2D, (|T]) further reduces to 



E 



E E E< 

i„+i=lm'G^„(Jf) JP'' 



(n) fem',rn+i 



V Ik I 

Z^ I "-ml 



(54) 
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exp{-(A+ Y. |k„|2^)(s-t)}P„(^)Q(^) n K^',Jt^^'dK 



n [i-e- 



^{|kj,|2^ + mp/3)(s-t) 



E 



K-',do,o,dk,)v^^,jo,o,dk,) 



^p i^p 



q i'9 



Here X^ denotes the summation over all incomplete Feynman diagrams JF' based on the set 
^„(J^) U {n + 1}. 



The proof of the lemma is complete with 






fl cr^' (ki, ■ ■ ■ , k„+i) := <fl(^{ki, ■■■, K) — ±r-^- 

' L-i I '^m' I 

m'eA„(j^) 



+ 00 



A+ E Ikm 



2,(3 



-(A+ ^ |k„|2^)(s-t) 



™eA(^') TT h g-(lkp|2'3 + lk,l2/')(s-t) 



Proof of Lemma §. 

We recall the pseudogenerator of the path corrector 



dSr 



E 



£x^")(t):=lim. 



^(n)(t + 5)|V_^^^ 






X^"Hi) 



5iO (5 

in the L^ sense. By a finite Taylor expansion we have 

(71),^ + 5 



^W(t + S)= e"x5H^, ^e(t)) + i?5(t), 



.£^ V ^2 



(55) 



with 



i?5(t)=£'^-^u(?)(-,a;,(t))5 + o(5) 



where o{6) is in the L^ sense. 
^From (p^) and (^5]) we have 

^"E x5^(^>^e(t))|V_,. 



+0O 

e" / e" 

t+s 



f2ro *+'5^„ 



ug-^)(s,^.W)|v„ 






(56) 
+ o{6). 



Thus, 



^Xi-^t) = X^-^t) - U^r'\t) + U(")(t) = Ceit; xS)) + U(")(t). (57) 

Here C^it] a;), t > denotes the pseudogenerator of the process e^x^^Wi x),t >0 calculated 
with respect to the filtration V_^ t , t > for a fixed a; G -R"'. The processes 



Mi")(t;a;) :=e"x^")(— 



x) — Ce{s; x)ds, for a fixed x G R"^ (58) 
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and 

t 

M^-\t) = xf\t) - Uxf\s)ds (59) 



are both martingales (cf. Section ^. 

The quadratic variation of the first martingale can be calculated with the help of Propo- 
sition H and it is given by 






where 



W^%{s,x)= (60) 



^2„-2 



Y.' j---j vtavt^cT'Pn-i{r)Q{r)Pn-i{r)Q{:F') 

a(im)|k/^ / KmK','^', 



m|2a+d-2(^2+ Y. |k^|2/3)(£2+ ^ \^J^P)\'^\' |k,|2 

^a„a;,'5(kg-k,/)t/kgt/kg/ \{ V^^^i^{^,x,dk^) Y[ V^:^^i>J^,x,dk'J. 

m^An{T),m^q m£A„(:F'),inj^q' 

The summation on the right hand side of the above expression extends over all multiindices 
of length n, i, i' whose first components are i and j correspondingly, all pairs of diagrams J-", 
T' G Qs made of vertices of 5 = {1, ■ ■ ■ , n}, all pairs of proper functions cr, cr' : 5 — i> {0, 1} 
and all pairs of g, g' belonging to A„(jF), A„(jF') respectively. We have denoted here k^ = 

As for the martingale given by (^9|), after elementary calculations we get 

t 
<M^f,M^:^>,= Jc{M^^M^:^}{s)ds (61) 



with 

£ {m^M^} is) = C {x^Ix^} is) - Cx^hs)x^]{s) - xS{s)Cx^{s). (62) 

Calculating precisely as in (|55|) we obtain that the first term on the right hand side of (|6^ ) 
equals 

Here £^jj(s;a3), s > denotes the pseudogenerator of the process e^^xT {'^■,^)Xj {-^^x)-, 
s > with X E W^ fixed. Taking into account (^) we obtain that 

£ {M^fM^-/} is) = Wi%{s, xM)- (63) 

In particular, the above calculation shows that M^^^ is a Brownian Motion whose covari- 
ance matrix equals (p6|). 
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4.1 Proof of eqs. ( |37|) - (^). 

We break the proof into two cases: a + 2/5 > 2 and a + 2(3 <2. 

The case of a + 2/? > 2. /^From a + /5 < 1 and a + 2/? > 2 it follows that /? > 1. 
For 2/? > 1, there exists a constant C{n,P), depending only on n and /3, such that, for any 
rrij G Aj{J^), 



mGAj(:r) 



m6Aj(JF) 



<c 



Ik^nl +e0 



«-mo I 



p2 _|_ IV 12/3 ■ 



So, 



"-1 IW I -I- p 



2/3 



(64) 



(65) 



J=l I "'J I 

for any rrij G Aj(jF). Let m-j := j if j is not the right endpoint of an edge. Otherwise, let 
nij be the closest free vertex to the left of the edge whose right endpoint is j. 
/^From (p5|) we have 



E 



in) 



xr(o) 



< 



(66) 



^,2„-2/,^ E /■■■//■•7 n 



i|2 J^<^Qs,cn ' 



i\kJ+E^){\k'J+E^] 



n 



(|kp| 



e2 + |k^|2/3)(^2+|k^|2/3) 

er){\k'j^\+et)a{\kp\)a{\k;\)dkpdk'p 



i+«„ 



VqdEir) 



{e^ + \K?'){^' + \V)^\K\K\) 



d+2a-2 



E n v^.„(o,rfk„)v;„(o,rfk:, 

where am > is the number of left endpoints between a free vertex m G A„(jF) and the 
next free vertex m' G A„(jF). Here and below the variables k'-s denote the second set of 
wavenumbers arising as a result of squaring of the right side of the expansion formula for 
-^C") (^^^. We also use the notation k„+p = k^, thinking that the second copy of the diagram 
is based on the shifted set {n + l,n + 2, ..., 2n}. 

Adopting the notation, kp = \kp\, subsequently and computing the expectation on the 
right side of (^) we get that 



K 



E 



J^ee. :f' 



K 

I n 

m£An{T) 



{kra + ei^){k'^ 



i. 



(e2 + kW) [e^ + {k'J^P) 



l+a„ 



(67) 



n 



0\f\/p fvfj ]Qjf£pQjfvQ 



n 



l.2a-l 

p n&E{r) 



{kp + e0){k'p + e^)dkpdk'p 
+ kf) (e^ + {k'pfP) {kpk'pY-^ 
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where T' is summed over QJy^AniTy) U {n + A„(jF)}. 
Note that 

c„ + 2 ^ oim = n, (68) 

where c„ denotes the cardinahty of the set A„(jF). 

For a + 2/5 > 2 the integrals with respect to the variables /Cp, p - the vertices of the left 
edges of JF', become unbounded when e | 0. Indeed as a result of integrating out the Dirac's 
delta functions we obtain terms whose singularity at are becoming, when £ J, of the same 
order of magnitude as (2+c.p+a,)^2/3-i)+2^-i >> tWss^- This fact gives rise to divergence 
of the corresponding improper integral at in the above regime of parameters a,/3. Since 
/? > 1, otherwise a + 2/? < 2 and we are not considering that case here this singularity can 
b Then changing variables by setting kp = kpe~'^ we obtain using ( p5D that 

K K 



E 



^:''''(0)V < CYY e"~^''^"''T''~"~'" I ... I TT {kp + e0){k'^ + eP)dkpdk'^ 






b pqeE{j^) \f + kp j [e + {k'^) pj {kpk'p) " 

(69) 



Ke '^ Ke ^ 

/ •/ n 



{kp + l)2+"p+"^dfc 



v 



r>,m^')[kr+if^'''^'''~kr-^' 

The first multiple integral appearing on the right hand side of (^9]) remains bounded when 
e I thanks to a + (3 < 1. The utmost right hand side factor however becomes then a 
product of improper integrals. The rate of decay of the integrand corresponding to index p 
in that integral is the same as for 

<< ~.^.o„. o , for L >> 1. 



r(2/3-l){2+ap+a,)+2a-l ^^ r.4/3+2a-3 ' """ '"P 
Kp tvp 

Since a + 2/3 > 2 the integrals involved are convergent. (^) follows then from the fact that 
c„ > 1. 

Proof of (|39|) . From Lemma ^ and (|25D we have that 



rn.gA„(Jf)U{n+l} 

where U^") = (t/i , ■ ■ ■ , t/^" ) and Iv?-^ | < 1- E|U„(0, 0)p can be estimated in the same 
way as x*'"'*- We outline the argument in the following. 

First (^) and the choice of rrij we get a bound, analogous to (1661) : 

e2(""i)E|U(")(0,0)|2< 
17 



i T 

n 



n 

meA„(:F)u{n+l} 



(|k™|+£^)(|k:„|+£^ 



(im + £^)(|k;i + £^)a(|kp|)a(|k;i)cikpdk; 



l+a„ 



Vq&E{r) 



e' + \K?'){^' + \V)^\K\K\) 



d+2a-2 



E 



n t/,„(o,cikj\/,„(o,rfk:,) 

,meA„(:F)u{n+l} 



Here we define a^ 



+1 



-1. Repeating the calculations leading to ( |69D and noting 



mGA„{:F)Un+l 



n-2. 



we obtain using 






n+en+2(en + l)(l-a-/3) 




Ke 1^ Ke T^ 






(71) 



(72) 



K 



_ft: 



n 






2a-l 



2/3 , i^2+ap+a,~2„__^ 



P<1 { K^ + 1 

with J^' G ^e(A„(J^) U{n + l}U{n + l + A„(J^) U {n + 1}}. 

The multiple integral appearing in the last factor on the right hand side of (173) becomes 
improper when we pass to the limit with e J, 0. However the integrals involving variables 
with both Op, aq nonnegative are convergent since a + 2(3 > 2. The only possible divergence 
comes therefore from the factors involving «„+! = — 1- The worst type of divergence could 
take place for the exponent On+i appearing twice in such a factor. This case corresponds to 
an occurrence of the edge (n + 1, 2n + 2). The end result is then the following asymptotic 

^2(n-i)g|u(n)(0,0)|2 < Cj2e 3 < Ce 3 . 

The assertion of the lemma follows for n > 2 [(3] . 

Proof of (|39|). We can use ( |63D to represent the quadratic variation of the martingale 



in question. After taking its expectation and using [|T0] we get the following estimate, valid 
for any p > 1 



E|MW(t)r<Ct^^(El^i;i(0,O)) 



>\P/2 



(73) 



i=l 



where VFg",; is defined by (|60|) . On the other hand using that formula one gets estimating in 
the precisely same way as it was done for (|37D and (^) that 

K K 



Eiyi3(o,o)^<c7.--EE E E /■■■/ 



(74) 



n 

m&A„{T) 






l+a„ 



(•^) ^' 

2n integrals 

k^J^' ( km' + e-p 

km' [kml + £ ) \ ^»T>' + ^ 



2a„ 



m' \ rn,' 



i,2o-l 11 



i,2o-l 



g ^ {kp + e'^){k'p + e'^)dkpdk'p 



pq(iE(T') 

with r' G ^c(A„(J^) U {n + A„(J^)} {m', m' + n}. 

Two situations arise: a^' = and a^' > 1- For a^' = 0, with an appropriate change of 
variables we rewrite the term in the above summations 

K K 



y...y_n 





{kp + ei3){k'p + ef>)dkpdk'p 
pqeE{j^) (^^ + ^p j (^£^^ + [kpr^j [kpk'pY"'-'^ 



K 



1,2/3 



_ 1 _ 1 



^m' I'^m' + ^ 



2^2 



-dkm.' 



n 

mgA„(:F') 






l+a„ 



n 

m(LE{r) 



01 /bp rhq jCirCpCirCn 



Ce" 



2(en-l)(l-a-/3) + Cn+n-2 



For Om' > 1, a similar calculation leads to 

K K 



£^"-^/---/ n 



{kp + ef^){k'p + e0)dkpdk'p 







fjEir) {^' + kf) {e^ + {k'pY^) {kpk'pf--^ 



_ 1 _ 1 

Ke '? Ke 7 






n 



{k^+i){k'^+i) 



l + Ori 



-p-p (5(A;p — kq)dkpdkq 2c„(i-a-/3)+c„+n-2 

11 p^^^l ^^-' ~ ^^ " 

pq&E{T') P 

In either situation, each term is of order 0{e"') for some 7 > 0, and Vn > 2. 

The special case of a + 2/5 = 2 can be similarly analyzed. This time, some integrals 
would diverge logarithmically, as e ^ 0, but would be controlled by factors of order 0{e'^) 
for some 7 > 0. 
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The case of a + 2/? < 2. We shall prove (|37|)-(|39D only for n = 1. This is sufficient 
for the proof of Theorem |l|. The general case requires evaluating multiple spectral integrals 
as in the previous case of a + 2/5 > 2, only much more easily because the corresponding 
improper integrals are convergent. 

We have 



E 



xi^Ht)r = E|x«(0)f<.^/^^,^d^< 



£2 + A;2/3)2 



Ce — 3 — 



dk 



[1 + A;2/3)2A;2"-i 



which vanishes, as e — ^ 0, for a + (3 < 1. Also, 



E 



t(i)/ 



U\^'(0,0) <CE|V(0,O)rE Vxr(0,O) 



,Wi 



and, consequently. 



limsupE Ul^^(0,0^ 



< +00 



as 



E 



Vxi'^(0,O) 



< 



i?<* 



|kpg(|k|) 

(A+|k|2/3)2 



2-a-2/3 



rfk<C(l + A^^) 



(75) 



(76) 



(77) 



is bounded, as A -^ 0, for a + 2j3 < 2. The inequality (^) holds for any p-th moment of 
Uj^ , with p > because of Gaussianity of the velocity field. 
We also have 

L 2 J- „ 



E 



JljfJ{±,x,{s))ds =2|rf.|E{uf;(^,a..(s))-U«(^,a;,(.i))}rf^i 







t s 



2jdsjE\v^^^i^,x,is,))-V^^^C^,xM))^ds,+ (78) 



s 

|rf.|rfsi|E|(v(|,^,(.2))-VUf.H^,a.,(s2)))-U(^2^(^,a..(5i))}rf32 



t s 



si 

Because of stationarity of U^ along the trajectory (see [|l^) the first term of (^) equals 

t s 



2 J ds jE{v^J^i^,O).vi'J{0,O)^ds, 







t s 



2|rfs|E|(v(^,0;0)-Vx^2^(^,0;0))-Uf2^(0,0)}rf.i. 



(79) 
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Here V(-, ■; 0) denotes the orthogonal projection of V(-, ■) onto L'^^q. A similar notation 
has been used for Vxe^- 

By incompressibility of the velocity 



r(i) 



(1)/ 



V\''it,x;0)=\it,x;0)-xrit,x;0). 
Using Cauchy-Schwarz inequality we estimate (|79|) by 

ft s 

JdsjE\Vxf^{'^,O;0)\'ds,\ . (80) 



1/2 



^0 



Since the fourth moment of U^ remains bounded as e | 0, (|80|) can be estimated by using 
(0) and (|73) as 

1/2 
C2 



.0 



^2/3 



X 



^2a-l(^2 + p/3y 



-dkds 



1) 



which tends to for a + 2/5 < 2, in view of the Dominant Convergence Theorem. 
The second term in (^) can be rewritten, after a simple change of variables, as 

t s s 

S2 

By the spectral representation of V we obtain 



VV^'^it^x-s) 



R^ Ri 






W(s, a;, rfk, dki) 



where 



W(s, a;, dk, dki) := k ■ Vi(s, x, c?ki)ki (g) V(s, a;, t/k) + k ■ V(s, a;, (iki)k (g) Vi(s, a;, d\s). 



Moreover by stationarity of the Lagrangian velocity experienced by the particle (see |10|) 
and incompressibility of velocity we have 

2 



rE 



jyv^^['-l^x,{s,)-'-^)ds, 



1 



rE 



S — S2 



rE 



l_,^p{_(W!!±^Sip(£z£2)- 



j VU^'.)(^,O;0)dsi 
} 



(83) 



R'i R'i 



(|k|2/3 + £2)(|k|2/3+|ki 



2/31 



■W(0,O,cik,dki 
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which equals 



1 — exp I - 



(|kp/? + |kip^)(s-S2) 
^2 



(|k|2/3 + e2)(|k|2/3+|k^|2/3) 



X 



1 _ p^T. ^ (|kT^+|k'iP''^)(s^S2) 



(|k'|2/3 + |k;|2/5)(|kf/5 + e2 



ij'* i?/' _Rd ^d 

E {W(0, 0, dk, (iki)W(0, 0, dk', dk'i)} . 
By incompressibihty of velocity again, the above expectation equals 

|ki|2trR(k)k ■ R(ki)k + |kpk • R(ki)R(k)ki) 6(k - k[)6{ki - k') + 



(|k|2trR(k)k ■ R(ki)k + |kipk • R(ki)R(k)ki) (5(k - k')5(ki - k^^ 
Thus we can bound the utmost left hand side of (|8^) by 



dkdkidk'dk[. 



C< 



K K 

II 





k^{k + ki 



1 — exp 



(fe2'3+fcf){s-S2) 

7S 



(F/3 + £2)2pa-1^2(/3+a)-l 



K K 



k\k + ki] 



(A;2/3 + ki 



1 — exp 



2/3 N 



dkdki + 



(84) 



(fc2/3+fcf )(s-S2) 
f2 







pa+/3-1^2a+/3-1^^2/3 + ^2) (^2/3 _^ ^2) 



(A;2/3 + A;?^) 



■dkdki 



with the constant C independent of e. For a + 2/5 < 2, (^) tends to zero, as e | 0, by the 
Dominant Convergence Theorem. Thus, the second term of ( [78| ) vanishes with en 

References 

[1] R. J. Adler: Geometry of Random Fields, Wiley, New York, 1981. 

[2] R. A. Carmona and J. P. Fouque: Diffusion approximation for the advection diffusion of 
a passive scalar by a space-time Gaussian velocity field, Seminar on Stochastic Analysis, 
Random Fields and Applications. (Edited by E. Bolthausen, M. Dozzi and F. Russo). 
Birkhauser, Basel, 37-50, 1994. 

[3] I. I. Gihman and A. V. Skorohod: Theory of Stochastic Processes I, Springer- Verlag, 
Berlin, 1974. 

[4] J. Glimm and A. Jaffe: Quantum Physics, Springer- Verlag, New York, 1981. 

[5] I. S. Helland: Central limit theorems for martingales with discrete or continuous time. 
Scan J. Statist. 9(1982), pp. 79-94. 



22 



[6] S. Janson Gaussian Hilbert Spaces, Cambridge Tracts in Math. 129, Cambridge Univ. 
Press, 1997. 

[7] H. Kesten, G. C. Papanicolaou: A limit theorem for turbulent diffusion, Comm. Math. 
Phys. 65 (1979), pp. 97-128. 

[8] R. Kubo: Stochastic Liouville equation, J. Math. Phys. 4 (1963), pp. 174-184. 

[9] H. J. Kushner: Weak Convergence Methods and Singularly Perturbed Stochastic Control 
and Filtering Problems, Birkhuser, 1990. 

[10] S. C. Port and C. Stone: Random measures and their application to motion in an 
incompressible fluid, J. Appl. Prob. 13 (1976), pp. 499-506. 

[11] M. Rosenblatt: Markov Processes. Structure and Asymptotic Behavior, Springer- Verlag, 
Berlin, Heidelberg, New York, 1971. 

[12] Yu. A. Rozanov. Stationary Random Processes, Holden-Day, 1967. 

[13] Yu. A. Rozanov. Markov Random Fields, Springer Verlag, 1980. 

[14] Ya. G. Sinai. On higher order spectral measures ofergodic stationary processes. Theory 
Probab. Appl. 8 (1963), pp. 429-436. 

[15] A. N. Schirayev. Some questions on spectral theory of higher moments. Theory Probab. 
Appl. (Russian) 5 (1960), pp. 295-313. 

[16] G.I. Taylor: Diffusions by continuous movements, Proc. London Math. Soc. Ser. 2, 20 
(1923), pp. 196-211. 



23 



